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f 要 : 依据 采 自 青海 省 杂 多 县 昂 赛 乡 的 树 轮 样本 ,建立 树 轮 宽度 差 值 年 表 ( RES)。 相 关 普 查 表明 ,澜沧江 源 区 


6 一 7 月 平均 的 12 个 月 时 间 尺 度 标准 化 降水 蔡 散 指数 (SPET) 与 RES 年 表 具 有 显著 的 正 相 关 。 利 用 RES 年 表 重建 


了 澜沧江 源 区 1639 一 2013 年 6 一 7 月 平均 的 SPEI 变化 序列 ,重建 方程 方差 解释 量 达 到 37. A96 , 且 经 过 交叉 检验 证 
明了 方程 的 稳定 性 。 重 建 序列 展示 了 澜沧江 源 区 干旱 的 长 期 变化 情况 ,发 现 1694 一 1710 .1753 一 1778 ,1830—1847 , 
1862—1872 .1892 一 1908 、1919 一 1929 、1978 一 1989 年 共 7 个 偏 湿 阶 段 和 1646—1673 , 1682—1693 , 1711—1731 , 
1735—1752 .1796 一 1809 .1817 一 1829 1848 一 1861 ,1873—1886 ,1949—1959 .1990 一 2000 年 10 个 偏 干 阶段 。 重 建 的 
SPEI 序列 存在 30 ~64 .19 ~21 4.2 ~4.3 2.6 -3.22.5 a 和 2.2 a 左右 的 显著 变化 周期 。 相 关 分 析 表 明 , 赤 道中 东 
太平 洋 海 温 对 澜沧江 源 区 干 湿 变 化 具有 显著 影响 。 此 外 ,澜沧江 源 区 干 湿 变 化 可 能 与 厄尔尼诺 - 南方 涛 动 (EN- 
SO) .太平 洋 年 代 际 振荡 (PDO) 、 太 阳 活 动 等 大 尺度 的 气候 驱动 因子 存在 一 定 的 联系 。 此 次 SPEI 重建 为 澜沧江 源 
区 提供 了 长 期 干旱 变化 历史 ,可 为 该 区 水 资源 管理 提供 科学 依据 。 
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三 江 源 地 处 青藏 高 原 腹地 ,是 长 江 、 黄 河 、 澜 沧 
江 的 源头 ,被 誉 为 “中 国 水 塔 "" 。 该 地 区 地 形 复 
杂 , 气 候 多 变 ,是 中 国 淡水 资源 的 重要 补给 区 ,其 生 
态 环境 和 气候 变化 不 仅 会 影响 当地 水 资源 的 开发 和 


料 具有 地 域 分 布 广 且 易于 采样 .定年 准确 .时 间 分 辩 
率 高 .与 气候 水 文 要 素 相关 度 高 等 优点 ,国内 许多 学 
者 利用 树木 年 轮 资料 研究 了 过 去 的 气候 变化 "”"。 
在 青藏 高 原 三 江 源 区 ,已 有 一 些 关 于 降水 .气温 、 径 


利用 ,而 且 对 我 国 乃 至 全 球 气候 变化 也 会 产生 重要 
的 影响 2。 澜沧江 发 源 于 青海 省 玉树 自治 州 ,位 
于 三 江 源 核心 区 ,干流 全 长 4 880 km ,流域 面积 约 为 
8.1 x 10? km? ,年 径流 量 4 633 x108 ms0G6) 。 流 域 水 
资源 主要 靠 雪 冰 融 水 和 大 气 降 水 补给 ,受气 候 变 化 
影响 较 大 。 近 年 来 ,全 球 气 候 变化 引起 了 流域 气温 、 
降水 ,径流 和 燕 发 等 气候 水 文 要 素 的 变化 ,使 流域 水 
资源 量 和 干支 流 流量 发 生变 化 ,导致 干旱 洪涝 等 极 
端 事件 频率 和 强度 增加 ,对 整个 流域 的 水 资源 环境 、 
自然 生态 系统 造成 严重 威胁 。 因 此 ,深入 研究 泣 
沧 江 源 区 干旱 历史 变化 特征 及 其 可 能 的 气候 影响 机 
理 , 预 测 未 来 干旱 变化 的 趋势 ,对 流域 生态 环境 、 经 
济 建设 以 及 水 资源 管理 具有 重要 的 科学 意义 。 

然而 ,由 于 该 区 域 仪 器 记录 的 气候 数据 年 代 较 
短 ,为 了 更 好 地 理解 区 域 气候 变化 的 长 期 历史 ,需要 
背 助 高 分 辩 率 的 古 气候 代用 指标 。 鉴 于 树木 年 轮 资 
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流 等 气候 水 文 要 素 的 重建 工作 ""“" , 但 目前 在 澜 沧 
江 源 区 基于 局 地 树 轮 资料 的 干旱 重建 还 不 多 见 。 
此 ,还 需要 进行 更 多 的 树 轮 气候 学 研究 ,以 获得 更 多 
关于 该 区 域 历史 干旱 变化 的 信息 。 
标准 化 降水 蔡 散 指数 SPEI (standardized pre- 
cipitation evapotranspiration index ) 是 在 标准 化 降水 
指数 SPI (standardized precipitation index ) 的 基础 
上 ,引入 潜在 燕 发 项 所 构建 "”。SPEI 是 一 种 较 新 
的 干旱 指数 , 既 具 有 多 时 间 斥 度 的 优点 ,又 考虑 了 气 
温 变 化 对 干旱 的 影响 。 它 不 仅 能 够 识别 干旱 事 
件 的 开始 和 结束 时 间 , 而 且 可 以 依据 干旱 的 强度 和 
持续 时 间 衡 量 其 严重 程度 ,目前 该 指数 在 全 球 范围 
已 被 广泛 应 用 于 干旱 分 析 和 评估 SPEL 的 数 
值 越 大 说 明 湿润 程度 越 高 ,数值 越 小 说 明 干 旱 程 度 
越 强 。Gou 等 “利用 树 轮 宽度 资料 重建 了 青藏 高 
原 东 北部 SPEI 过 去 千年 的 变化 历史 ,表明 SPEI 在 
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该 区 域 的 树 轮 学 研究 中 已 受到 关注 。 本 文 依据 采 自 
澜沧江 源 区 的 树木 年 轮 数 据 , 重 建 了 自 1639 年 以 来 
的 375 a 该 区 域 的 12 个 月 时 间 尺 度 6 一 7 月 平均 
SPEI 的 变化 历史 ,基于 重建 结果 ,分 析 区 域 干 旱 历 
史 变 化 特征 ,并 初步 分 析 区 域 干旱 变化 的 可 能 气候 
影响 机 理 。 以 期 为 政府 部 门 制定 适应 和 减缓 旱 洲 灾 
害 的 对 策 以 及 建立 合理 的 水 资源 可 持续 管理 方案 提 
供 新 的 基础 资料 。 


1 材料 与 方法 


1.1 研究 区 概况 和 资料 选取 

研究 区 地 处 青藏 高 原 东 北部 澜沧江 上 游 , 采 样 
点 位 于 青海 省 杂 多 县 昂 赛 乡 , 为 高 原 寒 温 性 气 
候 2 。 区 内 海拔 较 高 ,树木 较 少 受 人 类 活动 影 
响 55 。 在 遵循 树 轮 气 候 学 采样 点 选择 的 基本 原理 
基础 上 ,还 考虑 了 坡 向 与 空间 分 布 等 因子 的 影响 , 进 
行 合理 的 采样 布局 (图 1, 表 1)。 采 样 树种 为 大 果 
IRE (Sabina tibetica) , 主要 生长 在 海拔 3 992 ~ 
4 092 m, J BE Jy 30? ~40° 的 南 坡 上 ,采样 时 间 为 2013 
年 9 月 。 总 共 采 集 了 28 棵 树 ,得 到 56 AANE HR 
离 采 样 点 较 近 的 杂 多 (95°18'E,32°54'N) UL SE WE 
(96729'E,32?12' N) 两 个 气象 站 ,选取 的 时 间 范 围 均 
为 1961 一 2013 年 。 利 用 Mann-Kendell 方法 和 double- 
mass analysis 方法 对 气象 要 素 进行 均一 性 检验 '“-”， 


表 1 树 轮 采 样 点 概况 


Tab.1 Survey of the tree-ring sampling sites 


采样 点 ”北纬 ”东经 RES WE 坡 向 海拔 
名 称 (N) Œ) / 树 00 /m 


昂 赛 乡 32°42' 95?37' 56/28 30-40 PHJ 3992-4092 


36°N T T 


”采样 点 
。 气 象 站 
* SPEI 格 点 


35°N 上 


34°N F 


33°N 上 


32°N 


31°N : : : : 
93"E  94°E 95°E 96°E 97E 98E 
图 1 树 轮 采样 点 .气象 站 位 置 和 标准 化 降水 
蒸 散 指数 (SPEI) 格 点 分 布 


Fig.1 Locations of tree-ring sampling sites and meteorological 


stations and distribution of SPEI grid points 


结果 表明 两 站 的 气象 要 素 无 随机 突变 和 明显 的 不 均 
匀 分 布 。 由 图 2 可 见 , 杂 多 和 客 谦 站 月 平均 气温 和 
月 降水 量具 有 明显 的 季节 性 。 月 平均 气温 7 月 最 高 
1 月 最 低 , 月 降水 量 在 6 月 或 7 月 最 大 。 暖 季 (5 一 9 
月 ) 两 站 降水 量 分 别 占 全 年 降水 量 的 85. 9% 和 
88.5% 。 该 区 域 夏季 (6 一 8 月 ) 降水 量 最 大 是 雨 热 
同期 ,冬季 寒冷 干燥 。 由 于 多 个 气象 站 的 气候 要 素 
平均 值 比 单个 站 点 的 能 更 好 代表 区 域 气候 信息 , 因 
而 本 文采 用 杂 多 和 品 谦 两 站 的 算数 平均 序列 进行 后 
续 分 析 。 

本 文 使 用 的 SPEI 源 于 荷兰 皇家 气象 研究 所 的 
数据 共享 网 站 (http://climexp. knmi. nl/) 。SPEI 数 
据 集资 料 的 分 辩 率 为 0.5° x0.5°, 采 用 的 时 间 尺 度 
为 12 个 月 ,时 间 范 围 为 1961 一 2013 年 ,选取 距离 采 
样 点 最 近 的 SPEI 格 点 (32.75°N,95.75?°E) 进 行 分 
析 。 

1.2 年 表 研 制 

在 中 国 科 学 院 地 理科 学 与 资源 研究 所 树 轮 实 验 
室 完 成 了 树 轮 样本 的 量 测 ,并 严格 按照 目前 国际 通 
用 的 基本 程序 进行 了 样本 的 预 处 理 。 为 了 剔除 奇异 
点 过 多 和 与 主 序列 之 间 相 关 性 较 差 的 个 别 样本 , 利 
用 国际 树木 年 轮 库 的 COFECHA 程 序 对 交叉 定 

一 平均 气温 一 平均 最 高 气温 口 降 水 量 


[ (a 杂 多 站 


降水 量 /mm 
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图 2 1961—2013 年 杂 多 和 吉 谦 气象 站 月 降水 量 、 
平均 气温 和 月 平均 最 高 气温 的 年 内 分 布 


Fig.2 Mean monthly rainfall and temperature and maximum 


temperature at Zaduo and Nangqian meteorological 


stations during 1961 — 2013 


I) 


TK 


年 结果 进行 了 质量 控制 。 然 后 ,使 用 ARSTAN 4X 
EO 建立 年 表 。 通 过 对 原始 序列 进行 比较 分 析 ,并 
反复 试用 不 同 的 拟 合 方法 ,最 终 选 用 了 步 长 为 100 
的 样 条 函数 进行 树木 生长 趋势 的 拟 合 。 对 个 别 生 长 
趋势 不 符合 上 述 函 数 曲 线 的 序列 ,需要 调整 使 用 负 
站 数 函数 和 Hugershoff 曲线 "进行 分 别 拟 合 。 最 
后 ,将 拟 合 后 的 各 序列 合成 得 到 3 种 类 型 的 树 轮 宽 
度 年 表 : 标 准 化 年 表 ( STD) 、 差 值 年 表 ( RES) 和 自 回 
归 年 表 ( ARS)。 并 选择 含有 较 多 高 频 信息 的 RES 
年 表 进 行 后 续 分 析 。 表 2 给 出 了 采样 点 RES 年 表 
的 部 分 统计 特征 及 公共 区 间 分 析 结 果 , 由 表 2 可 见 ， 
RES 年 表 可 靠 时 段 为 1639 一 2013 年 , 共 375 a, 


表 2 采样 点 树 轮 差 值 年 表 ( RES) 统计 特征 及 公共 区 间 
(1800—2000 年 ) 分 析 结 果 
Tab.2 Statistical characteristics and common intervals 
(1800 —2000) of residual tree-ring chronology ( RES) 


昂 赛 乡 统计 量 
1474 一 2013 年 年 表 长 度 
1639/9 SSS > 0. 85 的 起 始 年 / 样 芯 数 
0.27 平均 敏感 度 
0.212 年 表 标 准 差 
0.40 所 有 序列 相关 系数 
0.40 树 间 相关 系数 
0.63 树 内 相关 系数 
21. 42 信 噪 比 
0. 96 样本 量 总 体 代表 性 (EPS) 
42.6% 第 一 主 成 分 方差 解释 量 


2 结果 与 分 析 


2.1 树 轮 宽 度 对 气候 因子 的 响应 
为 了 分 析 树 轮 年 表 所 包含 的 气候 信息 ,利用 距 
离 采 样 点 较 近 的 杂 多 、 守 说 两 个 气象 站 区 域 平均 的 


降水 和 气温 数据 ,将 RES 年 表 与 上 一 年 9 月 到 当年 
8 月 各 气象 要 素 序 列 进行 相关 分 析 。 结 果 表 明 ,RES 
年 表 与 当年 5 月 (P<0.05).6 月 (P<0.05) 和 8 H 
(P<0.01) 降水 呈 显 著 正 相关 ,同时 ,与 当年 5 月 
(P<0.01) 平 均 气 温和 5 月 (P<0.01).6 月 (P< 
0.05) fll 8 月 (P<0.05) 平 均 最 高 气温 呈 显 著 负 相 
关 ( 图 3)。 已 有 的 研究 ”表明 ,在 干旱 - 半 干 旱 
地 区 ,树木 生长 季 内 降水 量 减少 ,气温 的 升 高 会 导致 
树木 呼吸 和 蒸腾 作用 加 强 ,土壤 的 燕 散 作用 加 强 , 从 
而 限制 树木 的 生长 。 树 轮 宽度 指数 与 降水 正 相关 ， 
与 气温 负 相 关 的 气候 响应 模式 为 湿度 压力 型 生长 模 
HOD ,这 种 生长 模型 在 青藏 高 原 地 区 已 有 的 树 轮 研 
jg pais Sup 07 79 。 

考虑 到 SPEI 是 个 多 时 间 尺 度 的 干旱 指数 ,以 及 
RES 年 表 与 降水 和 气温 的 显著 相关 性 ,笔者 计算 了 
RES 年 表 与 各 月 时 间 尺 度 上 SPEI 的 相关 系数 , 结 
表明 ,RES 年 表 与 SPEI 指数 的 相关 系数 随时 间 尺 度 
的 不 同 而 变化 ,相关 系数 最 高 的 SPEI 时 间 尺 度 为 
12 个 月 ,说 明 研究 区 大 果 圆 柏 能 够 适应 短期 的 干 
i ,树木 的 生长 对 干旱 的 响应 主要 发 生 在 长 时 间 尺 
度 上 。 因 此 ,笔者 选取 了 12 个 月 时 间 尺 度 的 SPEI 
用 于 后 续 的 干旱 分 析 和 重建 。 如 图 4 所 示 ,在 12 个 
H RE E, RES 年 表 与 SPEI 数据 的 相关 系数 明显 高 
于 与 降水 和 气温 等 气象 因子 。 树 木 生长 与 上 一 年 9 
月 至 当年 8 月 连续 12 个 月 的 SPEI 都 存在 较 好 的 正 
相关 关系 ,其 中 以 与 当年 6 月 相关 系数 最 高 (r = 
0. 591 ,P <0. 001) ,与 当年 5 月 和 7 月 相关 系数 也 
达到 0.5 以 上 (图 4) 。 这 表明 水 分 是 限制 研究 区 大 
果 圆 柏 径 向 生长 的 主要 因素 , 而且 干 旱 指数 比 降水 
数据 能 更 好 的 代表 水 分 条 件 对 树木 生长 的 限制 。 对 
于 多 月 组 合 的 气象 要 素 ,RES 年 表 与 当年 6 一 7 HF 
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注 :p9 .p10 ,p11 和 pl2 分 别 代表 上 一 年 的 9 月 10 月 11 月 和 12 月 。 下 同 。 
图 3 采样 点 树 轮 宽度 指数 与 各 月 气候 要 素 之 间 的 相关 关系 


Fig.3 Correlation coefficients between tree-ring width index and climate factors in different months 
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图 4 树 轮 宽 度 差 值 年 表 (RES ) 与 SPEI 的 相关 关系 


Fig.4 Correlation coefficients between residual tree-ring chronology ( RES) and SPEI index 


均 的 SPEI 相关 性 最 好 ,相关 系数 达 0.611(P < 
0.001) 。 因 此 , 它 被 确定 为 用 于 重建 的 气候 要 素 。 
2.2 ”SPEI 的 重建 与 检验 
在 上 述 相关 分 析 基 础 上 ,采用 线性 回归 的 方法 
建立 了 澜沧江 源 区 12 个 月 时 间 尺 度 的 6 一 7 月 平均 
SPEI, ,与 RES 年 表 之 间 的 回归 方程 : 
SPEI, , =3.648RES - 3.228 


式 中 :SPEI。 ;代表 6 一 7 月 平均 澜沧江 源 区 SPEI 
RES 代表 树 轮 宽度 差 值 年 表 序 列 。 重 建 方程 的 相关 
系数 达到 0.611(n=53,P<0.01) ,重建 方程 可 以 解 
释 研究 区 6 一 7 月 SPEI 序列 方差 的 37.4% ,调整 自 
由 度 后 为 36.2% 。 为 验证 重建 方程 的 可 靠 性 ,笔者 
采用 逐一 别 除 法 对 方程 进行 检验 ,检验 的 各 项 统计 
量 列 于 表 3。 结 果 表 明 ,重建 序列 与 实际 值 符号 一 
致 的 年 份 有 35 a (通过 了 0.05 显著 性 水 平 检验 ) ， 
一 阶 差 符号 一 致 的 有 34 a (通过 了 0.01 显著 性 水 
平 检验 )。 此 外 ,重建 方程 的 乘积 平均 值 i 为 3. 03， 
误差 缩减 值 (RE) 为 0. 33 ,下 检验 值 为 30. 46 ,都 通 
过 了 0.01 的 显著 性 检验 ,上 述 统 计量 均 表明 了 重建 
方程 能 较 好 地 反映 低频 和 高 频 的 变化 旦 稳定 可 靠 。 
图 5 给 出 了 重建 的 6 一 7 月 SPEI 在 1961 一 2013 年 
与 实际 值 的 对 比 , 从 图 中 可 以 发 现 , 重 建 的 澜沧江 源 
区 1961 一 2013 年 6 一 7 月 SPEI 序 列 基 本 能 够 再 现 


lim 


A3 转换 方程 和 交叉 检验 的 统计 量 
Tab.3 Statistical parameters of the transfer 


equation and cross-validation 


一 阶 差 符号 mS 乘积 平均 值 a ”误差 缩 
检验 符号 检验 a) F WCRE) 
34(32*,34^) 35(35*,37") 3.03 30.46 0.33 


注 :a 表示 达到 0. 05 显著 性 所 需 同 号 数 ;b 表示 达到 0.01 显著 性 所 
需 的 同 号 数 。 


实际 的 SPEI 变化 特征 ,重建 曲线 所 反映 出 的 高 低频 
变化 特征 与 实际 曲线 都 具有 良好 的 一 致 性 。 

2.3 重建 序列 变化 特征 分 析 

通过 计算 得 出 ,重建 的 澜沧江 源 区 SPEL ;序列 
平均 值 为 0.372 ,标准 差 o 为 0.763。 故 定义 :高 于 
平均 值 + 1o 的 年 份 为 湿润 年 ,反之 , 低 于 平均 值 
-lo 的 年 份 为 干旱 年 。 在 过 去 的 375 a 中 ,研究 区 
的 湿润 年 有 58 a ,干旱 年 有 53 a ,分别 占 总 年 份 数 的 
15.596 和 14. 196 ,正常 年 有 264 a, 占 总 年 份 的 
70. 496 ;极端 湿润 年 份 ( > 平均 值 +2o) 有 5 a, 极 端 
干旱 年 份 ( < 平均 值 -2c) 有 12 a ( 表 4)。 与 青海 
玉树 地 区 过 去 500 a 的 早 涝 事件 进行 比较 可 见 ， 
重建 序列 结果 中 反映 的 极端 湿润 年 , 如 1703、1704、 


ar 一 重建 值 


一 实际 值 ， 


_3 pop pe i f , 
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图 5 1961—2013 年 SPEI 实际 值 与 重建 值 比较 
Fig.5 Comparison of the recorded and reconstructed values 


of SPEI during 1961 —2013 


表 4 1639—2013 年 澜沧江 源 区 极端 干 湿 年 
Tab.4 Extreme wet-dry years in Headwaters of the 
Lancang River during 1639 — 2013 
极端 湿润 年 1640 ,1699 ,1703 ,1704 ,1873 


1689 ,1700 ,1724 ,1727 ,1739 , 1749 ,1872 ,1910 , 
1942 ,1953 ,1995 ,1998 


极端 二 早年 
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1910,1942 .1953 .1995 年 也 与 旱 涝 级 别 中 的 少雨 年 
( 旱 涝 等 级 三 4) 相 一 致 ,进一步 说 明了 重建 序列 的 
可 靠 性 。 

为 研究 重建 序列 的 年 代 际 变化 特征 ,对 该 序列 
做 11 a 滑动 平均 处 理 后 (图 6) ,定义 连续 11 a 以 上 
高 于 均值 的 时 段 为 偏 湿 阶段 ,连续 11 a 以 上 低 于 均 
值 的 时 段 为 偏 干 阶段 ,发现 7 个 偏 湿 持续 阶段 和 10 
个 偏 干 持 续 阶段 ( 表 5)。 最 长 偏 湿 时 段 为 1753 一 
1778 年 , 共 持 续 26 a; 最 长 偏 干 时 段 为 1646 一 1673 
年 , 共 持 续 28 a, 

利用 多 窗 谱 分 析 方 法 “(MTM) 对 SPEI, ; 重 
建 序列 进行 周期 分 析 , 以 了 解 近 375 a 来 澜沧江 源 
es DA 性。 由 MTM 谱 估 计 和 其 相应 的 

声 临 界 谱 分 析 ( 图 7) 发 现 ,重建 的 6 一 7 月 平均 
pa i a 34 ~ 64.19 -21,4.2 -4.3,2. 6 ~ 
3.2.2.5 a 和 2.2 a 左右 的 显著 周期 变化 (达到 了 


RS 1639 一 2013 年 澜沧江 源 区 干 湿 时 段 
Tab.5 Wet - dry periods in Headwaters of the Lancang 
River during 1639 — 2013 


emm — C0. FMB A SpE 

1694—1710 年 0. 86 1646—1673 年 0.15 
1753—1778 年 0.56 1682—1693 年 0.12 
1830—1847 年 0.55 1711—1731 年 0.15 
1862—1872 年 0.55 1735—1752 年 0.20 
1892—1908 年 0.57 1796—1809 年 0.30 
1919—1929 年 0.55 1817—1829 年 0.28 
1978—1989 年 0.59 1848—1861 年 0.17 

1873—1886 年 0.20 

1949—1959 年 0.19 

1990—2000 年 0.28 


0 0.1 0.2 0.3 0.4 0.5 


图 7 SPEI 序列 的 MTM 谱 估计 及 其 相应 的 红 噪 声 临 界 谱 
Fig.7 MTM spectrum of SPEI series and its corresponding 


red noise critical spectrum 


9596 红 噪 声 置信 水 平 ) 。 其 中 34 ~64 a 和 19 ~21 a 
左右 的 周期 可 能 与 太平 洋 年 代 际 振荡 (PDO)" 及 
太阳 活动 有 关 *。4.2 ~4.3.2.6~3.2.2.5 a 和 
2.2 a 的 显著 周期 在 厄尔尼诺 — 南方 涛 动 (ENSO ) 
变 率 的 范围 内 "” ,表明 澜沧江 源 SPEI 变化 可 能 与 
ENSO 具有 一 定 的 联系 。 


3 ”区 域 比较 


为 了 检验 重建 SPEL ;的 空间 代表 性 ,分 别 对 重 
建 和 观测 的 SPEI。_， 序 列 与 同期 全 球 格 点 的 SPEI 数 
据 进行 空间 相关 分 析 。 结 果 显 示 ( 图 8a) ,重建 
SPEL ;与 同期 青藏 高 原 中 东部 广大 地 区 的 SPEI H. 
有 显著 正 相关 ,相关 系数 最 高 的 区 域 位 于 青海 南部 
高 原 澜 沧 江 源 区 (相关 系数 7 20.6, P <0. 001,n = 
53), 。 虽 然 与 观测 (图 8b) 相 比 ,重建 序列 的 相关 系 
数值 相对 小 一 些 , 但 两 个 空间 相关 场 分 布 非常 相似 。 
观测 和 重建 的 SPET ;与 我 国 青藏 高 原 大 部 分 区 域 
的 同期 SPEI 变化 均 有 显著 的 正 相 关 。 表 明 此 次 重 
建 的 SPEI, ;可 以 代表 较 大 区 域 范围 的 SPEL ， 
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注 :“ +” 表示 树 轮 采 样 点 位 置 ;SPEIe 7 代表 6 一 7 月 平均 澜沧江 源 区 SPEI。 深 ( 浅 ) 阴影 区 表示 通过 99% (9596 ) 的 显著 性 检验 。 下 同 。 


到 8 SPEI, ;重建 (a) 与 观测 (p) 序 列 的 空间 相关 性 


Fig.8 Spatial correlation between the reconstructed values (a) and observed series (b) of SPEI, , 


变化 。 

为 进一步 验证 重建 序列 的 可 靠 性 ,选取 了 附近 
地 区 已 有 的 研究 结果 进行 了 对 比 , 秦 宁 生 等 "重建 
的 青海 南部 高 原 春 季 (4 一 6 月 ) 湿润 指 数 序列 ,Xiao 
等 "重建 的 澜沧江 源 区 扎 曲 河流 域 年 流量 序列 , 秦 
宁 生 等 重建 的 长 江 源 区 通天 河流 域 年 径流 量 序 
列 ,时 兴 合 等 “重建 的 青海 杂 多 5 一 6 月 干燥 度 指 
数 序 列 。 经 过 计算 发 现 ( 图 9) ,本 文 重 建 的 SPEI, 
与 扎 曲 河 年 流量 .通天 河 径流 量 和 青 南 高 原 4 一 6 月 
湿润 指数 序列 呈 显 著 正 相关 且 都 达到 了 0.001 的 显 
著 性 水 平 , 相 关系 数 分 别 为 0.54 .0. 47 .0. 46 , 与 杂 
多 地 区 5 一 6 月 干燥 度 指 数 呈 显著 负 相 关 (r = 
-0.56,P <0.001)。 对 重建 序列 进行 11 a 滑动 平 
均 后 ,相关 系数 分 别 为 0.61.、0.38.、0.36 和 -0.42， 
均 达 到 0. 001 显著 性 水 平 ,表明 5 条 曲线 在 低频 变 
化 上 表现 出 较 好 的 一 致 性 。5 条 序列 在 18 世纪 00 
年 代 .19 世纪 30 年 代 到 40 年 代 初 .19 世纪 90 年 代 
末 到 20 世纪 00 年 代 初 为 一 致 的 湿润 期 ,以 及 17 世 
纪 80 年 代 末 到 90 年 代 初 .18 世纪 40 年代、19 世纪 
50 年 代 末 到 60 年 代 初 为 一 致 干旱 期 。 这 说 明 澜 沧 
江 源 区 干 湿 变化 具有 较 强 的 一 致 性 。 

此 次 重建 的 SPEI。 ;所 展示 的 年 代 际 干 湿 阶 段 
变化 ,也 可 以 在 青藏 高 原 地 区 重建 的 干 湿 序 列 中 得 
到 一 些 印 证 。 例 如 ,本 序列 中 存在 的 1735 一 1752 
年 .1873 一 1886 年 .1949 一 1959 年 几 个 偏 干旱 时 段 
与 青藏 高 原 东北 部 5 一 6 月 PDSI 序列 中 显示 的 
干旱 持续 阶段 较为 一 致 。 重 建 的 SPEI 序列 中 19 t 
纪 20 年 代 的 偏 干旱 持续 阶段 与 多 个 基于 树 轮 重建 
的 青藏 高 原 地 区 干 湿 序 列 -较为 一 致 。 这 种 


干 湿 阶 段 的 一 致 性 ,不 仅 证 明了 此 次 重建 的 SPEI, ; 
的 可 靠 性 ,还 表明 了 青藏 高 原 澜沧江 源 区 的 干 湿 变 
化 对 大 尺度 气候 变化 的 响应 在 低频 变化 上 具有 一 定 
区 域 代表 性 。 另 外 ,将 重建 的 SPEI。 ;序列 与 Cook 
等 "重建 的 亚洲 地 区 夏季 帕 尔 默 干旱 指数 ( Mon- 
soon Asia drought atlas, MADA ) 格 点 数据 进行 对 比 ， 


1.00 (a) 本 文 重建 的 澜沧江 源 6 一 7 月 SPEI 指 数 
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图 9 11 a 滑动 平均 后 的 几 个 重建 序列 对 比 


Fig.9 Comparison of 11 year running averages between our 


reconstructed SPEI series and those in previous studies 


TOS 区 


计算 重建 序列 与 MADA 中 澜沧江 源 区 采样 点 距离 
较 近 的 格 点 重建 数据 在 1639 一 2005 年 期 间 的 相关 
系数 ;结果 显示 ,重建 的 SPEL ;序列 与 附近 的 4 个 
MADA #4 (93. 75?E,33. 75? N93. 75?E,31. 25?N; 
96. 25^E,33. 75? N96. 25^? E,31. 25? N) JE TA 4H 2C Z% 
数 都 达到 了 0. 001 显著 性 水 平 ,分 别 为 0. 488, 
0.389 0.391 和 0.221 ,由 此 进一步 证 明了 此 次 重建 
序列 的 可 靠 性 。 


4 与 全 球 海 温 的 关系 


海 温 的 变化 可 以 通过 影响 大 尺度 大 气 环流 系统 
异常 ,进而 影响 区 域 的 降水 , Fc e b 79 。 
为 了 进一步 研究 澜沧江 源 区 干旱 变化 的 气候 机 理 ， 
基于 上 述 对 重建 序列 的 周期 分 析 结 果 , 将 重建 的 澜 
沧 江 源 区 12 个 月 时 间 尺 度 的 6 一 7 月 SPEI 序列 与 
前 一 年 7 月 至 当年 6 月 平均 的 全 球 海 表 温 度 场 NC- 
DCv3b ERSST “进行 了 空间 相关 分 析 ,初步 分 析 重 
建 的 澜沧江 源 区 干 湿 变 化 与 海 温 的 关系 。 

相关 分 析 结 果 如 图 10 所 示 , 在 公共 时 段 
1961—2013 年 期 间 ,重建 的 SPEI。 ,序列 与 赤道 中 
东 太 平 洋 (10。$ ~ 109N,120? ~ 180°W) 海 温 时 显著 
负 相 关 ; 与 热带 西 太 平 洋 (10° ~ 20? N, 130? ~ 
150°E) 海 温 呈 显著 正 相 关 。 男 外 ,在 南半球 太平 洋 
西南 部 也 具有 显著 的 正 相 关 区 域 。 当 赤道 中 东 太 平 
洋 海域 海 温 偏 高 ( 低 ) 时 ,澜沧江 源 区 SPEL -7 MAR 
(高 ) ,澜沧江 源 区 气候 偏 干 旱 ( 湿 润 ) ; 当 热 带 西 太 
平 洋 区 域 和 南半球 太平 洋 西 南部 海 温 偏 高 ( 低 ) 时 ， 
澜沧江 源 区 SPEL -M r R) ,澜沧江 源 区 气候 偏 

60°N 


20°E 60°E 100°E 140°E 


180° 
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湿润 (干旱 )。 表 明太 平 洋 (尤其 是 赤道 中 东 太 平 洋 ) 
海 温 对 澜沧江 源 区 干 湿 变 化 具有 明显 的 影响 ,ENSO 
对 澜沧江 源 区 气候 干 湿 变 化 的 影响 具有 季节 延续 性 。 
这 与 吕 爱 锋 等 “对 三 江 源 地 区 径流 变化 对 ENSO 信 
号 的 响应 关系 结果 以 及 Xiao 等 "9 对 澜沧江 上 游 扎 
曲 河流 域 流量 与 太平 洋 海 温 的 相关 结果 一 致 。 

为 进一步 探讨 区 域 干 湿 变化 与 大 尺度 环流 系统 
的 联系 ,将 重建 的 SPEI。 ;与 上 一 年 7 月 至 当年 6 月 
的 NINO 3. 4 ( http://www. esrl. noaa. gov/psd/data/ 
climateindices/list) 指数 序列 进行 公共 时 段 内 
(1951—2013 4E) 的 相关 分 析 ,发 现 二 者 存在 显著 负 
相关 关系 (r= -0.33,P <0.01,n=63)。 进 一 步 表 
明 澜 沧 江 源 区 干旱 变化 与 ENSO 有 密切 联系 。 已 有 
的 基于 观测 资料 的 研究 -指出 ,ENSO 对 东亚 气 
候 的 影响 是 通过 一 个 正 的 热 动 力 海 气 反馈 机 制 而 实 
现 的 。 当 海 温 处 于 ElNino 位 相 ,赤道 中 东 太 平 洋 海 
温 升 高 会 引起 大 气 Rossby 波 响 应 ,有 利于 西北 太平 
洋 海 温 偏 冷 ,导致 在 菲律宾 地 区 出 现 一 个 反 气 旋 式 
环流 异常 ,使 东亚 夏季 风 偏 弱 , 季 风 输 送 到 我 国 北方 
的 水 汽 明 显 减 弱 , 从 而 引起 澜沧江 源 区 气候 偏 干 旱 。 
而 当 海 温 处 于 LaNifia 位 相 , 情 况 则 相反 。 

此 外 ,笔者 将 重建 序列 进行 11 a 滑动 平均 处 理 
后 ,发现 其 与 PDO 指数 (http://www. atmos. washing- 
ton. edu/ ~ mantua/abst. PDO. html) 显著 相关 ,相关 
系数 达 -0.44 (P «0.01) ,暗示 了 PDO 和 青藏 高 原 
地 区 气候 变化 的 联系 ,但 PDO 变化 周期 较 长 ,而 现 
有 的 气象 观测 资料 年 限 较 短 ,其 对 研究 区 干旱 变化 
的 具体 影响 机 制 仍 有 待 于 我 们 结合 古 气候 和 现代 气 
修 的 研究 方法 进行 深入 研究 。 


140°W 100°W 60°W 20°W 


图 10 1961—2013 年 澜沧江 源 区 SPEL, ,重建 序列 与 前 一 年 7 月 至 当年 6 月 海 温 场 相关 系数 分 布 
Fig. 10 Correlation coefficients between the reconstructed SPEI。， and the sea surface temperature field from July of last year 


to June of that year during 1961 —2013 
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5 结论 


(1) 在 青海 省 杂 多 县 昂 赛 乡 28 棵 大 果 圆 柏 上 
采集 了 56 个 样 芯 ,得 到 一 组 树 轮 宽度 样本 ,经 过 样 
本 预 处 理 交叉 定年 .标准 化 等 步骤 ,最 终 建立 了 
RES 年 表 。 和 气候 响应 分 析 表 明 , RES 年 表 与 澜沧江 


源 区 干旱 数据 SPEI 具有 良好 的 相关 性 , 且 具 有 显著 
的 树木 生理 学 意义 。 综 合 考虑 RES 年 表 与 气温 、 降 


JK SPEI 的 相关 性 ,选取 6 一 7 月 平均 12 个 月 尺度 
SPEI 用 于 干旱 重建 ,建立 转换 方程 并 单 点 重建 了 澜 
沧 江 源 区 近 375 a 来 干旱 历史 变化 ,经 过 交叉 检验 ， 
证 明了 此 次 重建 序列 的 稳定 性 和 可 靠 性 。 

(2) 重建 的 澜沧江 源 区 SPEI, ;展现 出 了 干旱 
的 长 期 变化 情况 , 发现 1694—1710 、1753 一 1778 、 
1830—1847 ,1862—1872 ,1892—1908 , 1919—1929 , 
1978—1989 年 共 7 个 偏 湿 持续 阶段 和 1646—1673 、 
1682—1693 .1711 一 1731 .1735 一 1752 ,1796—1809 、 
1817—1829 ,1848—1861 ,1873—1886 , 1949—1959 , 
1990—2000 4E 10 个 偏 干 持续 阶段 。1640 、1699 、 
1703 .1704 年 和 1873 年 为 极端 湿润 年 , 共 5 a;1689、 
1700 , 1724 , 1727 , 1739 , 1749 , 1872 , 1910 、1942 、 
1953 1995 年 和 1998 年 为 极端 干旱 年 , 共 12 a, 3E 
建 序 列 与 青藏 高 原 中 东部 干 湿 变 化 有 较 好 的 一 致 
性 , 且 在 年 代 际 时 间 尺 度 上 与 青藏 高 原 区 域 的 多 个 
重建 干 湿 序列 体现 的 干 湿 阶 段 较为 吻合 。 说 明 澜 沧 
江 源 区 的 干 湿 变 化 在 低频 上 对 大 尺度 气候 变化 的 响 
应 ,具有 一 定 的 区 域 代表 性 。 

(3) 周期 分 析 结 果 显 示 ,澜沧江 源 区 干旱 变化 
具有 30 ~64 19 ~21.4.2 ~4.3.2.6~3.2.2.5 a 和 
2.2 a 左右 的 显著 变化 周期 。 周 期 变化 也 体现 了 
PDO 、 太 阳 活 动 ENSO 变化 可 能 与 澜沧江 源 区 干 湿 
变化 具有 一 定 的 联系 。 相 关 分 析 表 明 , 赤 道中 东 太 
平 洋 海 温 对 澜沧江 源 区 干 湿 变 化 具有 显著 影响 。 
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Drought Reconstruction Based on Tree-Ring Width in Headwaters 
of the Lancang River 


HUANG Xiao-mei', XIAO Ding-mu’, QIN Ning-sheng? 
(1. Chengdu Institute of Plateau Meteorology , CMA/Sichuan Provincial Key Laboratory for Heavy Rain and Drought-Flood 
Disasters in Plateau and Basin, Chengdu 610072 , Sichuan, China; 
2. Sichuan Provincial Meteorological Observatory, Chengdu 610072 , Sichuan , China; 
3. Sichuan Provincial Climate Center, Chengdu 610072 , Sichuan, China) 


Abstract;  Tree-ring width residual chronology ( RES) was developed from the tree-ring samples collected in An- 
gsai Township, Zaduo County, Qinghai Province. The results of general survey showed that there was a high corre- 
lation between the tree-ring width chronology and the SPEI (June-July) at 12-month time scale. The June-July 
SPEI in the headwaters of the Lancang River was reconstructed from AD 1639 to 2013, and it was used to explain 
the 37. 496. of the total SPEI variance during their common period. The cross-validation revealed that the recon- 
structed equation was stable. The results showed that the relatively wet periods occurred during 1694 — 1710, 
1753 - 1778, 1830 - 1847, 1862 —1872, 1892 - 1908 , 1919 —1929 and 1978 — 1989, and the relatively dry pe- 
riods occurred during 1646 —1673, 1682 - 1693 , 1711 - 1731 , 1735 - 1752, 1796 —1809, 1817 —1829, 1848 - 
1861, 1873 —1886, 1949 —1959 and 1990 - 2000, respectively. The reconstructed SPEI series revealed that there 
were the significant variation periods of 30 264, 19 -21, 4.2 -4.3, 2.6 - 3.2, 2.5 and 2. 2 years. Correlation 
analysis showed that the sea surface temperature in the middle-east equator Pacific Ocean affected significantly the 
dry-wet changes in the headwaters of the Lancang River. Our reconstruction was associated with the large-scale at- 
mospheric-oceanic variability, such as the El Nifio-Southern Oscillation ( ENSO ), Pacific Decadal Oscillation 
(PDO) and solar activity. This SPEI reconstruction could provide a long-term drought history in the headwaters of 
the Lancang River and serve as the scientific basis for water resources management in this semiarid region. 


Key words:  Lancang Hiver; tree-ring width; SPEI reconstruction 


